Efficient readout of micromechanical resonator arrays in ambient conditions 
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We present a method for efficient spectral readout of mechanical resonator arrays in dissipative 
environments. Magnetomotive drive and detection is used to drive double clamped resonators in the 
nonlinear regime. Resonators with almost identical resonance frequencies can be tracked individually 
by sweeping the drive power. Measurements are performed at room temperature and atmospheric 
pressure. These conditions enable application in high throughput resonant sensor arrays. 
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Micro- and nanometer scale mechanical resonators are 
widely considered as mass and force sensors. Ad- or des- 
orption of molecules by the resonator are detected as 
slight changes in resonance frequency. Several methods 
exist to measure the resonance frequency of a single res- 
onator [l|, H, H, Small arrays of static mechanical 
sensors be scanned optically 0] or the reflection of the 
array as a whole may be analyzed @. Spectral readout 
of small-scale resonant arrays has been reported in vac- 
uum (3]. In this letter we present an efficient method 
to readout resonators with closely spaced resonance fre- 
quencies in dissipative environments. We show that the 
number of resonant sensors which can be operated in a 
given bandwidth is not limited by the Q factor. Using 
a magnetomotive setup we can discriminate mechanical 
resonators with almost identical frequencies in ambient 
conditions. 

To track two individual resonators with quality fac- 
tor Q in a frequency spectrum, their resonance frequen- 
cies fi should be separated by at least A/ = fi/Q. Simul- 
taneous operation of n resonators requires a minimum 
bandwidth on the order of A/ = /i(l -I- 1/Q)". Here /i 
is the resonance frequency of the slowest resonator. In 
vacuum, the Q-factor of a mechanical resonator can be of 
the order 10^ and readout of large arrays of resonators 
should be feasible. At atmospheric pressure however, Q- 
factors usually drop below 10^ because of viscous drag 
Q. This limits the number of resonators within a prac- 
tical bandwidth to a few hundred at most. Large-scale 
arrays therefore require either multiple detectors or sep- 
arated drive circuits, and the complexity of such systems 
rapidly increases with the number of resonators. 

To address this problem, we propose to readout ar- 
rays of resonators with closely spaced center frequen- 
cies by sweeping the drive power from low to high val- 
ues at constant frequency. Once the nonlinear regime is 
accessed, individual resonators are marked by instanta- 
neous transitions in the phase and amplitude response. 
When each nonlinear resonator is used as a sensing ele- 
ment 9] , the locations of these characteristic amplitudes 



are affected by ad- or desorption induced mass or stress 
change to the concerning resonator. We conduct exper- 
iments on arrays of double clamped resonators at room 
temperature and atmospheric pressure. We experimen- 
tally show that in a 2-D array of 6 resonators with closely 
spaced resonance frequencies, individual resonators can 
be discriminated using a single drive and detector unit. 
Calculations confirm the experimental findings. 

Arrays of double clamped mechanical resonators are 
fabricated out of 100 nm thick low-pressure chemical va- 
por deposited (LPCVD) silicon nitride by electron beam 
lithography. All resonators have the same dimensions 
200/im X 15/im. Slight variations in center frequencies 
of the resonators are expected as a result of variations 
in residual stress [l3|. Figure 1(a) shows one of the res- 
onators before metallization. A layer of chromium and 
gold is evaporated on top to enable magnetomotive drive 
and detection. The resonators share a support on one 
side, while supports at the other side are separated. This 
allows the measurement of both the individual and the 
collective responses. 

Figure 1 (b) shows the measurement setup schemat- 
ically. The RF voltage from a network analyzer is con- 
verted to a current which drives the resonators (two in 
series depicted). A magnetic field of 1.9 T is generated by 
a Halbach array pLlj constructed from NdFeB permanent 
magnets. The resulting voltage is amplified and mea- 
sured by the network analyzer. This measurement setup 
can be used at room temperature and atmospheric pres- 
sure, in contrast to earlier magnetomotive experiments in 
vacuum and/or at cryogenic temperatures. As a result 
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FIG. 1: (a) Scanning electron micrograph of a single res- 
onator, (b) Magnetomotive measurement setup. The current 
source is formed by the amplifier U/I and the resistor network. 
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FIG. 2: Magnetomotive measurements in ambient environ- 
ment, (a) Individual linear frequency response of two res- 
onators with almost identical resonance frequencies, (b) The 
same resonators driven nonlinear by sweeping drive power at 
215kHz. (c,d) Collective response of the resonators connected 
in series in linear (c) and nonlinear (d) regimes. (e,f) show 
the difference between the calculated sum of the individual 
responses and the measured collective response in the linear 
(e) and the nonlinear (f) regime. 



of viscous drag, the Q-factors in our experiments are a 
factor of 10 to 100 lower than in those experiments. 

Figure 2(a) shows the hnear response of two res- 
onators with almost identical center frequencies. The 
voltage drop due to the resistance of the resonator, 
bond wires etc. has been subtracted. The reso- 
nance frequencies and Q-factors are determined by fit- 
ting Lorentzian functions (shown in the figure) to the 
measured data. The resonance frequencies of these de- 
vices were 207.72kHz and 209.05kHz, and the Q-factor is 
39 for both resonators. The corresponding bandwidth is 
5.4kHz, whereas the difference in resonance frequency is 
1.3kHz. Now the resonators are connected in series and 
the collective response is measured. Figure 2(b) shows 
the result for the same drive conditions as in (a). Dis- 
crimination between the two resonators is impossible as 
the resonator bandwidth is more than four times the dif- 
ference in center frequency. We now individually drive 
the same resonators at increasing power at a constant 
frequency of 215kHz. Fig. 2(d) shows the result: a steep 
transition in the resonator amphtude marks the char- 
acteristic drive amplitude at this frequency. When the 
collective response is measured at strong driving, as in 
(d), the two resonators are easily distinguished. 

In absence of interaction, the response of multiple 
resonators is just the sum of the individual responses. To 
verify this, the difference between the sum of the individ- 
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FIG. 3: Measurements on a 2-D array of 6 resonators, (a) 
Resonator response while sweeping the drive current at / = 
270kHz. The inset shows the arrangement of the resonators 
and their linear resonance frequencies, (b) and (c) show am- 
plitude and phase response of the array, driven at different 
frequencies fd- 



ual responses of (a) and (d) and the collective response in 
(b) and (e) is plotted in Fig. 2 (c) and (f). No amphtude 
dependency is found, which indicates that the coupling is 
weak. Note that a high driving power results in a notably 
lower noise level in the nonlinear driven system p^ . 

We now turn to a 2-D array of 6 resonators, for 
which the arrangement and the measured linear reso- 
nance frequencies are depicted in the inset of Fig. 3(a). 
In this experiment, the bandwidth of the individual res- 
onators was 5 ±0.1 kHz, which corresponds to a Q-factor 
of approximately 50. 

Figure 3(a) shows the nonlinear amplitude response 
at fd = 270kHz; clearly 6 resonators can be distin- 
guished. In Fig. 3(b) and (c), sweeps at fixed drive fre- 
quencies ranging from 250kHz to 285kHz are plotted. 
Six lines marking characteristic jumps in resonator am- 
plitude and phase remain equidistant. Indeed, over the 
nonlinear regime, the total bandwidth covering the 6 res- 
onators changes by less than 60Hz in this measurement. 
When sweeping frequency from low to high values at con- 
stant drive current (not shown), we found strong coupling 
between the resonators which causes multiple resonators 
to collapse simultaneously at high driving power. Strong 
coupling makes the system useless as an array of inde- 
pendent sensors. 

As for the differences in the center frequencies of 
the resonators, we note that in LPCVD silicon nitride 
residual stress variations of the order of 10% across a 
2' wafer are not unusual [l0| . Using this number, the 
residual stress gradient would be on the order of 100 
Pa/^m. Given the distance between the resonators (at 
least 300/xm) the difference in residual stress is on the or- 
der of 10^ Pa which results in frequency differences on the 
order of lOOHz. Similar frequency differences can be ob- 
tained by slightly varying the geometry of the resonators 
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generated by the string-like linear resonator at resonance 
equals: 
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FIG. 4; Calculated amplitude responses, (a) Normalized am- 
plitude response of a single resonator driven at fixed frequency 
as a function of the normalized drive amplitude fd/fo- (b) 
Cross section at fd/fo = 1-2. (c) and (d) show the amplitude 
response for ten resonators. 



or the thickness of the deposited gold layer. 

To corroborate the experimental results, we have 
calculated the characteristic driving amplitude of a non- 
linear resonator at a given driving frequency by solving 
the equation of motion: 

My + Cy + {Kt + Kt,)v + K^y'^ = BLh cos{iut), (1) 

where M is the resonator effective mass, y the beam dis- 
placement C the damping constant, and Kt and Kb rep- 
resent linear stiffness terms due to residual tension and 
bending rigidity respectively. The nonlinear term K^ is 
a result of the elongation of the beam, L is the beam 
length, B the magnetic field strength and the driving 
current. 

For the beams in Fig. 2, the average residual tension 
To = 74/iN is found by comparing the calculated stress- 
free resonance frequency to the measured resonance fre- 
quency [l^. We calculated that Kt/Kh « 85, which 
indicates that the bending rigidity of the resonators is 
insignificant and the linear resonance frequency is merely 
determined by M and Kt- The magnetomotive voltage 
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With the experimental values B = 1.9 T, and Q = 39, 
and L = 200 /xm, the transduction at resonance equals 
Vemf/IcI - 0.8 V/A. 

Instead of the frequency-dependent resonator am- 
plitude, we now calculate the resonator amplitude while 
varying the drive amplitude at constant frequency fd to 
mimic the experimental situation. Thus, out of two sta- 
ble states, the resonator always vibrates at the lowest 
amplitude. Figure 4(a) shows the calculated response 
for a resonator similar to ones used in the experiments. 
The frequency axis is normalized to the linear resonance 
frequency /o, and traces for the normalized drive fre- 
quencies fd/fo from 0.5 to 2.0 are combined in (a). Red 
regions correspond to a high amplitude. The trace for 
fd/fo — 1-2 is shown in Fig. 4(b). Sweeping parallel to 
the horizontal axis, locations with multi- valued resonator 
amplitudes are accessible at large drive amplitudes, when 
fd > fo- Figure 4(c) shows the simulated response of an 
array of 10 resonators described by Eq. [2 where a slight 
diff'erence in Kt results in different characteristic ampli- 
tudes for each resonator. We assume a difference in linear 
resonance frequency for adjacent resonators equal to one 
tenth of the resonator bandwidth. Panel (d) shows a 
cross section of the amplitude response at fd/fo = 1-2, 
where /q is the linear resonance frequency of the slowest 
resonator. The phase response (not shown) also displays 
steep transitions similar to the amplitude, whose location 
can be determined with high accuracy hj. The calcula- 
tions thus reproduce the experimental findings in detail. 

In conclusion, we have demonstrated an efficient 
way to operate and readout mechanical resonator ar- 
rays in dissipative environments. The center frequen- 
cies of the linear resonators can be spaced very closely 
as each resonator is marked by its characteristic drive 
amplitude at a fixed frequency, instead of the linear res- 
onance frequency. As the exponential relation between 
required bandwidth and the number of resonators is cir- 
cumvented, this technique may enable the construction 
of high-throughput resonant sensor arrays. 
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